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The standard quantum search algorithm lacks a feature, enjoyed by many classical algorithms, 
of having a fixed-point, i.e. a monotonic convergence towards the solution. Here we present two 
variations of the quantum search algorithm, which get around this limitation. The first replaces 
selective inversions in the algorithm by selective phase shifts of ^ . The second controls the selective 
inversion operations using two ancilla qubits, and irreversible measurement operations on the ancilla 
qubits drive the starting state towards the target state. Using q oracle queries, these variations 
reduce the probability of finding a non-target state from e to e^'+^, which is asymptotically optimal. 
Similar ideas can lead to robust quantum algorithms, and provide conceptually new schemes for 
error correction. 



I. INTRODUCTION 

The quantum search algorithm is like baking a souffle . . . you have to stop at just the 
right time or else it gets burnt. Q/ 

The quantum search algorithm 3l consists of an iterative sequence of selective inversion and 
diffusion type operations. Each iteration results in a fixed rotation (which is a function of the 
initial error probability) of the quantum state in a two-dimensional Hilbert space formed by 
the source and the target states. The iterative procedure keeps the rotation going forever at a 
uniform rate. If we choose the right number of iteration steps, we stop very close to the target 
state, else we keep on going round and round in the two-dimensional Hilbert space. To perform 
optimally, therefore, we need to know the precise number of iteration steps, which depends 
upon the initial error probability or equivalently the fraction of target states in the database. 
When we have this information, the quantum algorithm leads to a square-root speed up over 
the corresponding classical algorithm for many applications, including unsorted database search. 
When this information is missing, we can estimate the required number of iteration steps using 
various "amplitude estimation" algorithms 's','?], but that necessitates an overhead of additional 
queries. When the total number of queries is large, the additional queries do not cost much, but 
when the total number of queries is small, the overhead can be unacceptably large. 

In this article, we address the problem of finding an optimal quantum search algorithm in 
situations where, (i) we do not know the initial error probability (perhaps only its distribution 
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or a bound is known) , and (ii) the expected number of queries is small (so that every additional 
query is a substantial overhead) . Such situations occur in pattern recognition and image analysis 
problems (where each query is a significant cost), and in problems of error correction and 
associative memory recall (where the initial error probability is small but unknown). We look 
for variations of the quantum search algorithm that ensure amplitude enhancement, and also 
outperform classical search algorithms. 

The strategy is familiar from classical computation, i.e. construct a quantum algorithm that 
"converges" towards the target state. That is, however, impossible to do by iterating a non- 
trivial unitary transformation; eigenvalues of unitary transformations are of magnitude 1, so 
the best that can be achieved by iterating them is a "limit cycle" and not a "fixed point". As 
described above, this is indeed what happens in case of the quantum search algorithm. To obtain 
an algorithm that converges towards a fixed point, some new ingredient is needed, and several 
possibilities come to mind: 

(a) Some property of the current computational state offers an estimate of the distance to the 
target state. This estimate can be used as a parameter to control the extent of the next iterative 
transformation, e.g. the Newton-Raphson method to find the zeroes of a function. If such an 
estimate is not available, we have to look for a different method. 

(b) Suitably designed but distinct operations are performed at successive iterations. Such a 
method can converge towards a fixed point even using unitary transformations. We illustrate 
this method using the Phase- 7r/3 search algorithm described in detail in Section II. 

(c) Irreversible damping is introduced in the algorithm without explicit use of any property of 
the target state. With the right type of irreversibility (i.e. when all eigenvalues of the fixed 
iterative transformation are less than 1 in magnitude), the algorithm converges, e.g. the Gauss- 
Seidel method for solving a set of linear algebraic equations. Within the framework of quantum 
computation, such an irreversibility can be introduced by projective measurement operations, 
and we provide such an algorithm in Section III Q . 

In a broader context, fixed point algorithms possess two attractive features by construction, 
which are not commonplace in generic quantum algorithms: 

(1) The initial state is guaranteed to evolve towards the target state, even when the algorithm 
is not run to its full completion. 

(2) Any errors due to imperfect transformations in earlier iterations are wiped out by the sub- 
sequent iterations, as long as the state remains in the problem defining space. 

These are powerful motivations to explore quantum fixed point algorithms, even if they have 
other limitations. 

Let us consider an unsorted database in which a fraction / of items are marked, but we 
don't have precise knowledge of /. We run a particular algorithm which has to return a single 
item from the database. If the returned item is a marked one, the algorithm has succeeded, 
otherwise it is in error. Without applying any algorithm, if we pick an item at random, then the 
probability of error is e = 1 — /. The goal of the algorithm is to minimize the error probability, 
using the smallest number of oracle queries. If / is sufficiently small, then we can use the 
optimal quantum search algorithm to obtain a marked item using 0{\/^/J) queries. A few more 
queries to estimate / or to fine-tune the algorithm is not a problem, because overall we gain a 
quadratic speed-up compared to the classical case requiring exhaustive search. But when / is 
large, the number of oracle queries is small, and the quantum search algorithm doesn't provide 
much advantage — in particular the rotation can overshoot the target state. In such a situation, 
a simple classical algorithm (select a random item and use a query to check if it is a marked 
one) may outperform the quantum search algorithm. The same considerations apply to the more 
general amplitude amplification algorithms 0. There the initial quantum state is a unitary 
operator U applied to a given source state |s), and the probability of getting a target state after 
measuring this initial state, |?7tsp, is analogous to /. The probability of error, which has to be 
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FIG. 1: In the standard quantum search algorithm (left), the state vector may overshoot the target 
state. In the algorithm of Section II (right), the state vector always moves towards the target state. 



minimized using the smallest number of queries, is the probabihty of getting a non-target state 
after measurement, e = 1 — |C/tsp. 

In Section II, we show that by replacing the selective phase inversions in quantum search by 
suitable phase shifts we can get an algorithm that always produces amplitude amplification. As 
shown in Fig.l, the change in phase shift eliminates overshooting and the state vector always 
moves closer to the target state. By recursively applying the single iteration derived for any 
unitary operator U, we develop an algorithm with multiple applications of U that converges 
monotonically to the target state. The optimal value of the phase shift turns out to be 7r/3, and 
we refer to this algorithm as the "Phase- 7r/3 search" Remarkably, the algorithm converges 
using reversible unitary transformations and without ever estimating the distance of the current 
state from the target state. Explicitly, the initial error probability e changes to e"^ in one 
iteration. Recursive application of the basic iteration n times requires qt = 3qi-i -I- 1 (go = 0) 
oracle queries at the i"^ level, and makes the error probability e"^ . Thus e = is the fixed point 
of the algorithm, and the error probability decreases as e^'^+i as a function of the number of 
queries q. 

Section III describes a different implementation of the fixed point quantum search, where 
irreversible measurement operations direct the current state towards the target state @. In 
it the inversion and diffusion operations are controlled in a special way by two ancilla qubits 
and their measurement. The same transformation is repeated at every iteration, but since the 
transformation is made non- unitary by measurement, the quantum state is able to monotonically 
converge towards the target state. A major advantage of this implementation is that the e^"^"*"^ 
convergence is achieved for all positive integer values of q, compared to the restricted numbers 
q-n = (3" — l)/2 for the Phase- 7r/3 search. 

Note that the best classical algorithm can only decrease the error probability as e"^^^ (not e"^, 
since the last iteration does not need a query). Thus the fixed point quantum search improves 
the convergence rate by a factor of 2. In Section IV, we illustrate this difference with an explicit 
example. 



II. THE PHASE- 7r/3 SEARCH 



Consider applying the transformation 

0{e,(t)) = URlU^RfU , 



(1) 
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to the state |s), where i?^ and i?f are selective phase shift operators for the source and the 
target state respectively. Note that if we were to choose the phase shifts as 9 ~ (f) — tt, we would 
get one iteration of the amplitude amplification algorithm 01 • 

In the following, we study the particular case = = 7r/3. We show that when the U 
operation drives the state vector from \s) to \t) with a probability (1 — e), i.e. \Uts\^ = (1 ~ e)j 
then 0(7r/3,7r/3) drives the state vector from the source state to the same target state with 
a probability of (1 — e^). The deviation from \t) hence falls from e to e"^. The striking aspect 
of this result is that it holds for any kind of deviation from \t). Unlike the standard quantum 
search algorithm, which would overshoot the target state when e is small (Fig.l), the Phase- tt/B 
search always moves towards the target state. This property is useful in developing quantum 
algorithms that are robust to variations in the problem parameters. 

Connections to error correction are already evident in the previous paragraph. Let us say 
that we want to drive a system from |s) to \t) state/subspace, using the transformation U, with 



driving probability \Ut, 



(1 — e). Then the composite transformation 0{tt/3, tt/B) will reduce 



the error from e to e"^. This technique is applicable whenever the transformations U, U\ Rs and 
Rt can be implemented. That will be the case when the errors are either systematic or slowly 
varying, and the transformation U can be inverted with exactly the same error, e.g. when the 
errors are due to environmental degradation of some component. The technique would not apply 
to errors that arise as a result of sudden disturbances from the environment. Quantum error 
correction is carried out at the single qubit level in many implementations, where individual 
errors are corrected by identifying the corresponding error syndrome. With the machinery of 
this paper, errors can be corrected without ever needing to identify the error syndrome. 



A. Analysis 

Now let us analyze the effect of the transformation 0(7r/3, tt/S) when it is applied to the state 
Is). We want to show that 



(2) 



Applying the operations, U , U\ Rt^^ and J7, one by one, to the state |s), wc find that 



0(7r/3,^/3)|s) - U\s) (^e"/3 + \Uts? (e"/^ - l)'^ + \t)Uts (e"/' - l) 



(3) 



The deviation of this superposition from the target state \t) is given by the absolute square of 
its overlap with the non-target states. This probability is 



(ti|0(^/3,7r/3)|5) 



|0(7r/3,^/3)t(l-|t)(i|)0(7r/3,V3)|s) 



t/3 



- 1 



Substituting \Uts\ = (1 ~ e), the above quantity becomes 



(4) 



(5) 



We now use this result recursively to construct the quantum algorithm for searching in presence 
of uncertainty. 
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B. Recursion 

A few years after its invention, the quantum search algorithm 0, was generahzed to a 
much larger class of applications known as the amplitude amplification algorithms . In 
these algorithms, the amplitude to go from the state |s) to the state \t) by applying a unitary 
operation U, can be amplified by repeating the sequence of operations, Q — IgU^ItU, where Is 
and /( denote selective inversions of the \s) and \t) states respectively. Note that the amplitude 
amplification transformation with four queries is 

U {isU^ItU) {isU^ItU) {isU^ItU) {LUUtU) . (6) 

When the operation sequence IgU'^ItU is repeated ry times, the amplitude in the U^t) state 
becomes approximately 2rjUts provided "qUts ^ 1- The quantum search algorithm is a particular 
case of amplitude amplification, with U being the Walsh-Hadamard transformation and |s) being 
the |0) state with all qubits in the |0) state. The selective inversions enable the amplitudes 
produced by successive iterations to add up in phase. Consequently, the amount of amplification 
grows linearly with the number of repetitions of Q, and the probability of obtaining \t) goes up 
quadratically. 

Just like the amplitude amplification transformation, it is possible to recurse the transfor- 
mation URgU^RtU to obtain larger rotations of the state vector in a carefully defined two 
dimensional Hilbert space. The basic idea is to define transformations Ui by the recursion 

U^+i = U^RsUlRtU^ , Uo = U . (7) 

Unlike amplitude amplification, there is no simple structure in the recursion, and it is not easy 
to write down the actual operation sequence for Ui with large i, without working out the full 
recursion for all integers less than i. Let us illustrate this for 1/2- 

Uo = U , C/i = UoRsU^RtUo = URsU^RtU , (8) 

U2 = UiRsUlRtUi = (URsU^RtU) Rs {URsU^RtUf Rt (URsU^RtU) 
= {URsU^RtU) Rs (u^'rIurIu''^ Rt {URsU^'RtU) 

^ U {RsU'^RtU) {RsU^R\u^ {RlU^RtU) [RsU"^ RtU) . (9) 

This is not as simple as the corresponding transformation for amplitude amplification, Eq.®. 

Now it is straightforward to show that if \Uts\'^ = 1 ~ then |(J7i)tsp = 1 — e^^ . The recursion 
obeyed by the number of queries is, g^+i = 3qi + 1, qo — 0, with the solution qi = (3* — l)/2. 
As a function of the number of queries, therefore, |(C/i)ts. P = 1 — e^^i+i^ 'pjjg error probability 
thus falls as e^'^"'"^ after q queries, which is an improvement over the classical algorithm where 
the error probability is e"^"*"^ after q queries (more details in Section IV). 

C. Fixed Point of the Algorithm 

First, note that both the amplitude amplification algorithm, Eq.®, and the phase shift 
algorithm, Eq.(|5J), perform selective operations for the state and so from an information 
theoretic point of view there is no confiict in having fixed points. However, unitarity would be 
violated if there was accumulation of the target state due to repetition of a fixed transformation. 
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In the amplitude amplification algorithm, exactly the same transformation is repeated and so 
unitarity does not permit any fixed point. Although the phase shift algorithm is very similar 
to the amplitude amplification algorithm, the transformation applied at each step is slightly 
different from the others, due to the presence of the four operations Rg, Rt, Rl, RI- It therefore 
gets around the unitarity condition that prevents the amplitude amplification algorithm from 
having a fixed point. 

The e^'^+i performance of the Phase- 7r/3 algorithm has been shown to be asymptotically (i.e. 
in the limit e — + 0) optimal 1^. It can be shown that application of the general operator 0{9, 

to \s) changes its component along \t±) by the scale factor exp ( ^ (0 — 0) ) — 4 sin | sin ^ | C/ts | 

In the asymptotic limit, \Uts\ 1, and the scale factor is minimized by the choice 9 ~ (p ~ I 



III. MEASUREMENT BASED SEARCH ALGORITHM 



To obtain a fixed point quantum search, we have to find an algorithm that successively 
decreases the probability of finding a non-target state. Let us say that the initial state of a 
quantum register is U\s) = sinajt) -I- cosaji^) — an unspecified superposition of the target \t) 
and the non-target \ti_) states. We can construct a fixed point search algorithm by repetitively 
measuring the state — whenever the oracle query identifies the measured state as the target state, 
we stop the algorithm, otherwise we use the diffusion operation to rotate the non-target state 
towards the target state. To implement this idea, we attach to the quantum register an ancilla bit 
in the initial state |0). The oracle query flips the ancilla when the register is in the target state. 
Now if we measure the ancilla, outcome 1 tells us that we are done, and the register will give us 
the target state. Outcome tells us that the register is in a superposition \tj_) of the non-target 
states, and its probability is the initial error probability e — cos^ a. To decrease this probability, 
we apply the diffusion operator UlgU'^ to the register, conditioned on the measurement outcome 
being 0. That reflects |t±)|0) about U\s)\0) to give the state sin2a|t)|0) + cos 2a|t_L) |0). The 
error probability has thus decreased by the factor cos^ 2a. Iterating the sequence of oracle query 
and diffusion operations n times, the error probability is reduced to cos^ a cos^" 2a — e{2e— 1)^". 
This convergence is better than the Phase- 7r/3 search when e > 1/3, but worse when e < 1/3. 

Our goal is to find an algorithm which gives optimal convergence for all values of e, without 
knowing any bounds that e may obey. We know that the rotation provided by a single iteration 
of the quantum search algorithm overshoots the target state when e is small. In particular, when 
e is less than 1/2, the overshooting is so large that the new state has a smaller overlap with 
the target state than the initial state. That gives an inkling for constructing a better directed 
quantum search — somehow set a lower bound of 1/2 for e. The natural lower bound for e is 0, 
but we can easily make it 1/2 by diluting the database with extraneous states. Explicitly, we 
add to the quantum register a second ancilla in the state |-|-) = (|0) -I- |l))/\/2, and make the 
oracle query conditional on this ancilla. This logic produces the following g-iteration algorithm: 

• Attach two ancilla qubits in the |0) state to the source state |s), i.e. |s) |0)|s)|0). (In 
what follows, we refer to the former ancilla as ancilla-1 and the latter ancilla as ancilla-2.) 

• Apply _ff ® [/ ® / to this extended register to prepare the initial state of the algorithm, 

• Iterate the following two steps, q times: 

Step 1: If ancilla-1 is in the state |1), then perform an oracle query that fiips ancilla-2 
when the register is in target state. 

Step 2: Measure ancilla-2. If the outcome is 1, the register is certainly in the target state, 
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Iterate q times 



ancilla-1 |0) 



register \s) 



ancilla-2 |0) 




1: exit 

FIG. 2: Quantum circuit for our g-iteration algorithm. 



SO exit the iteration loop. If the outcome is 0, then apply the joint diffusion operator 

{H (g) U)Iqs{H (E) Uy to the joint state of ancilla-1 and the register. 

• After exiting or completing the iteration loop, stop the algorithm and measure the register. 
The quantum circuit for this algorithm is shown in Fig. 2. 



A. Analysis 



Let us analyze the algorithm step by step. The initial state is 



IV'^) = (if OC/(g) J)|0)|s)|0) 



|0) 



V2 



sin a|t) + cos ajfj.) ) |0) 



(10) 



The initial error probability is e = cos^ a:, and the initial success probability is / = 1 — e = sin^ a. 
We work in the joint search space of ancilla-1 and the register, denoted by the subscript j, where 
only the state \tj) = \l)\t) acts as the target state and all the states |0)|t), |0)|tj_), |l)|f_L) act 
as non-target states. Let \t'j) denote the superposition of all the non-target states. In the joint 
search space, the initial state is 



sin a.\ , 1,,,,, 



where the unit vector \tj) is 



, , , , ^ / sm a , „, , , cos a , 
l*^)=^( ^|0)|i) + ^|0)|U) 



V2 



V2 



V2 



(11) 



(12) 



2 1 1/2 ^ ^2/{l+e). For later reference, 

note that the error probability after measuring the joint state \tj), i.e. the probability of finding 



and the normalization constant N is (cos^ ^ + ^ sin^ a) 



the register in the non-target state \t±), is 

I (0| {ti. I*;.) p + I m^\t'j)\^ = cos^ a = Nh. 



(13) 
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Step 1 of the algorithm, using an oracle query, flips ancilla-2 when the joint register state is 
\tj). In step 2, we measure ancilla-2. If the outcome is I then we stop the algorithm, because 
the register is in the target state. The probability of getting I is sin^a/2 = //2; we have 
effectively put an upper bound of 1/2 on the success probability using ancilla-1. If the outcome 
is 0, which has probabihty then the joint state is \t'j). In this case, we apply the joint 

diffusion operation using the joint source state \sj) = |0)|s). The joint diffusion operation is a 
reflection about (7J (8) J7)|0)|s) = Uj\sj), in the two-dimensional Hilbert space spanned by \tj) 
and \t'j) as shown in Fig. 3. The state Uj\sj) makes an angle a^, defined by sin^ aj = sin^ a/2, 
with the state \t'j). So reflecting \tj) about Uj\sj) gives us a state that makes an angle 2a j with 
1^^). Its component in -direction is cos2aj = 1 — 2sin^Q:j — 1 — sin^ a = e. Thus the joint 
diffusion operation produces the final state 

IV/) = f/jlt;) = VT^lW + e\t'^) . (14) 

After measuring |'0/), the probability of getting \tj) is e^, so the total probability of getting \t'j) 
after one iteration is N~'^e'^ . Iterating the algorithm will keep on decreasing this probability by 
a factor of at each iteration, and after q iterations it will become N~'^e^'^. So the net error 
probability after q iterations is (cf. Ea. (|I3|l ) 

= N^eiN-h^"^) = e^^+i , (15) 

which agrees exactly with the corresponding result for the Phase- 7r/3 search. 
The above analysis allows us to also deduce the following features: 

(1) e = 1 can be made a fixed point of the algorithm, instead of e = 0, by effectively interchanging 
the roles of I i) and|tj_). This is achieved by flipping ancilla-2, only when the joint state is |l)|i_L). 
Then the probability of finding the register in the state \t), after q iterations, becomes (1 — e)^'?+^. 
Note that the same behavior can be obtained in the Phase- tt/B search, by replacing either R^^^ 
with i?j or -Rs^^ with Rs ^ ^"^ , This e = 1 fixed point can be useful in situations where certain 
target states are to be avoided, e.g. in collision problems. 

(2) When we use fixed point quantum search to locate the target state in a database, the initial 
error probability is e = 1 — /. The number of oracle queries required to reduce this probability 
to o(l) obeys 

(1 - = 0(1) e-(2«+i)^ = o(l) . (16) 

Thus we need q = 0{l/f) oracle queries to find the target state reliably. This scaling of fixed 
point quantum search is clearly inferior to the 0{\/\/J) scaling of the quantum search algo- 
rithm ^2]. Still, as discussed in the introduction, fixed point quantum search can be useful in 
situations where / is unknown and large. 

(3) A practical criterion for stopping the iterative algorithm would be that the error probability 
becomes smaller than some predetermined threshold eth- Provided we have an upper bound 
e < Cup < 1, we can guarantee convergence by choosing < eth- In the above algorithm 
2g -I- 1 can take all odd positive integer values, which is less restrictive than the values of the 
form 3" allowed by the Phase- 7r/3 search. 

(4) Depending on the outcome of the ancilla-2 measurements, the algorithm can exit the iteration 
loop before completing it. So the average number of queries is always less than the maximum 
number of iterations q. This is better than the Phase- 7r/3 search, where all the queries required 
by the number of iterations must be executed. 

(5) It is possible to stick to unitary operations throughout the algorithm and postpone measure- 
ment till the very end. In such a scenario, the unmeasured ancilla-2 has to control the diffusion 
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Measure ancilla 
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Uj\sj) 



Diffusion 



Uj\sj) 



Measure register 



FIG. 3: Step-by-step quantum state evolution in our g-iteration algorithm. The double arrows relate 
the equivalent states in the original search space (left) and the joint search space (right). 



operation and all the subsequent iterations (i.e. they are executed only when ancilla-2 is in the 

|0) state), and it cannot be reused in the iteration loop. We need a separate ancilla-2 for every 
oracle query, to ensure that once the states \tj) and \t'j) are separated by an oracle query in an 
iteration, they are not superposed again by subsequent iterations. The whole set of q ancilla-2 
can be measured after the iteration loop, in a sequence corresponding to the iteration number, 
to determine the target state. In this version, the unitary transformation is different for each 
iteration (because each iteration involves a different ancilla-2 and different controls), and the 
unitary iterations converge to a fixed point. 
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IV. QUANTUM SEARCHING AMIDST UNCERTAINTY 

The original quantum search algorithm is known to be the best possible algorithm for exhaus- 
tive searching [3. llClj| . therefore no algorithm will be able to improve its performance. However, 
for applications other than exhaustive searching for a single item, we demonstrate that suitably 
modified algorithms may indeed provide better performance. 

Consider the situation where a large fraction of the states are marked, but the precise fraction 
of marked states is not known. The goal is to find a single marked state with as high a probability 
as possible in a single query. For concreteness, say the marked fraction / is uniformly distributed 
between 75% and 100%. In that shown below, the probability of error for the fixed point 

search algorithm is approximately one fifth of that for the best (possible) classical algorithm, 
and approximately one fifteenth of that for the best (known) quantum algorithm. 

Classical: The best classical algorithm is to select a random state and see if it is the target 
state (one query). If yes, return this state; if not, pick another random state and return that. 
The probability of error is equal to that of not getting a single target state in two random 
picks, i.e. = (1 — /)^, which lies in the interval [0,1/16]. The overall error probability is 
approximately 2.1%. 

Quantum: The best quantum search based algorithm available in the literature, for this 
problem, is by Younes et al. 0. Using q queries, it finds the target state with a probability 

/ (^^4^^^ + TH^), where cos/3 = e (see Eq.(59) of Ref. [ill). When q = 1, the success 

probability becomes / (l + 4(1 — /)^), which lies in the interval [15/16, 1]. The overall error 
probability is approximately 5.7%, which is worse than that for the best classical algorithm. 

Fixed point search: If we use the transformation 0(7r/3, 7r/3), with U being the Walsh- 
Hadamard transform and the source state \s) being the |0) state, then the probability of being in 
a non-target state after one query becomes = (1 — /)^, which lies in the interval (0, 1/64). The 
overall error probability is approximately 0.4%, far superior to the two cases described above. 

How the quantum algorithms perform is graphically illustrated in Fig. 4. For the particular 
distribution of / considered above. Fig. 4(a) shows that the graph of the fixed point search 
algorithm lies entirely above the graph of Ref. The difference between the two becomes 

even more dramatic if we consider multiple query algorithms, as in Fig. 4(b). As mentioned 
earlier, for the type of problems discussed in this section, the fixed point algorithm is the best 
possible quantum algorithm [^. 

V. CONCLUSION 

The variant of quantum search discussed in this article supplements the original search algo- 
rithm, by providing a scheme that permits a fixed point and hence moves towards the target state 
in a directed way. The fixed point property makes the quantum search robust, and naturally 
leads to schemes for correction of systematic errors 12]. 

We have presented two implementations of the fixed point quantum search, both optimal 
in the worst-case scenario. At first sight, requirements of monotonic convergence and unitary 
quantum evolution appear to be in confiict. But the conflict has been avoided by a recursive 
approach in the Phase- 7r/3 search algorithm and by irreversible projections in the measurement 
based search algorithm. Between the two, the measurement based search algorithm is better 
behaved than the Phase- 7r/3 search, because (a) it allows all positive integer number of oracle 
queries instead of a restricted set of integers, and (b) its average-case requirement of number of 
oracle queries is smaller. 

Deeper insights in to the mechanisms of these algorithms are still desirable. For instance. 
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FIG. 4: Comparison of the performance of the Phase- 7r/3 search algorithm with that of Ref.|ll||. when 
the fraction of marked states, /, varies between and 1: (a) after one query, (b) after 13 queries. 



what is so special about phase shift of 7r/3, and how does the change of phase shift (from tt 
to 7r/3) convert the amplitude amplification algorithm to something totally different? Can one 
improve the algorithms in the non-asymptotic region by making the phase shifts (or rotation 
angles) a function of the query sequence, in a manner reminiscent of local adiabatic quantum 
algorithms? Such questions are definitely worth exploring. 

The potential applications of these algorithms are in problems with small but unknown initial 
error probability and limited number of oracle queries, where they assure monotonic power-law 
convergence. For such problems, the fixed point algorithms are better than the best classical 
search algorithms by a factor of 2. Applications to pattern recognition and associative memory 
recall are under investigation. 



[1] G. Brassard, Searching a quantum phone book, Science 275 (1997) 627. 

[2] L.K. Grover, Quantum Mechanics helps in searching for a needle in a haystack, Phys. Rev. Lett. 
78 (1997) 325. 

[3] G. Brassard, P. Hoyer, M. Mosca and A. Tapp, Quantum amplitude amplification and estimation, 
AMS Contemporary Mathematics Series Vol. 305, eds. S.J. Lomonaco and H.E. Brandt, (AMS, 
Providence, 2002), p. 53, e-print quant-ph/0005055 



12 



Quantum algorithms with Gxcd point: The case of database search 



[4] L.K. Grover, Quantum computers can search rapidly by using almost any transformation, Phys. 

Rev. Lett. 80 (1998) 4329. 
[5] L.K. Grover, Fixed-point quantum search,, Phys. Rev. Lett. 95 (2005) 150501, e-print 

quant-ph/0503205 

[6] T. Tulsi, L.K. Grover and A. Patel, A new algorithm for fixed point quantum search. Quant. Inform. 

and Comput. (2006), to appear, e-print quant-ph/0503205 
[7] L.K. Grover, From Schrddinger's equation to the quantum search algorithm, Pramana 56 (2001) 

333, e-print quant-ph/0109116 
[8] S. Chakraborty, J. Radhakrishnan and N. Raghunathan, Bounds for error reduction with few quan- 
tum queries, Proceedings of APPROX- RANDOM 2005, LNCS 3624 (Springer, Berhn, 2005), p. 245. 
[9] C.H. Bennett, E. Bernstein, G. Brassard and U. Vazirani, Strengths and weaknesses of quantum 
computing, SIAM Journal on Computing 26 (1997) 1510. 
[10] C. Zalka, Graver's quantum searching is optimal, Phys. Rev. A60 (1999) 2746. 
[11] A. Younes, J. Rowe, J. MiUer, Quantum search algorithm with more reliable behavior using partial 

diffusion, quant-ph/0312022 
[12] B.W. Reichardt and L.K. Grover, Quantum error correction of systematic er rors using a quantum 
search framework, Phys. Rev. A72 (2005) 042326, e-print |quant-ph/05062'42] 



